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Abstract. Recently developed methods for PT-symmetric models can be applied to 
quantum-mechanical matrix and vector models. In matrix models, the calculation of all 
singlet wave functions can be reduced to the solution a one-dimensional PT-symmetric 
model. The large-N limit of a wide class of matrix models exists, and properties of 
the lowest- lying singlet state can be computed using WKB. For models with cubic and 
quartic interactions, the ground state energy appears to show rapid convergence to the 
large-N limit. For the special case of a quartic model, we find explicitly an isospectral 
Hermitian matrix model. The Hermitian form for a vector model with 0(N) symmetry 
can also be found, and shows many unusual features. The effective potential obtained 
in the largc-N limit of the Hermitian form is shown to be identical to the form obtained 
from the original PT-symmetric model using familiar constraint field methods. The 
analogous constraint field prescription in four dimensions suggests that PT-symmetric 
scalar field theories are asymptotically free. 
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1. Introduction 

Since the initial discovery of PT symmetry pj, there has been considerable progress 
in expanding both the number of PT-symmetric models and our knowledge of their 
properties [21 13]. However, models with continuous internal symmetry groups have not 
been extensively developed. Of course, the field theories relevant to modern particle 
physics have continuous symmetries, and it is natural to seek PT-symmetric models 
with similar continuous symmetries. Here we review recent progress we have made in 
the the construction and analysis of PT-symmetric models of scalars with 0{N) or U (N) 
symmetry [HE]. Most of the results will deal with quantum mechanical models, usefully 
regarded as one- dimensional field theories. Of particular interest is the construction of 
the large- iV limit, as this has proven to be a very powerful theoretical tool in the analysis 
of many different field theories. 

Hermitian matrix models appear in many contexts in modern theoretical physics, 
with applications ranging from condensed matter physics to string theory. Interest in 
the large- iV limit of matrix models was strongly motivated by work on the large- A'c limit 
of QCD [6], but interest today is much wider. For example, Hermitian matrix quantum 
mechanics leads to a construction of two-dimensional quantum gravity coupled to c = 1 
matter [7|. It is surprising that the construction of PT-symmetric matrix models is 
somewhat easier than the the construction of models with a vector symmetry. The 
matrix techniques pioneered in [8] for Hermitian matrix quantum mechanics can be 
extended to PT-symmetric matrix quantum mechanics. In these models, the contours 
of functional integration over the matrix eigenvalues are extended into the complex 
plane. The large- limit can then be taken in PT-symmetric matrix theories just as in 
the Hermitian case. Quantities of interest such as the scaled ground state energy and 
scaled moments can be calculated using WKB methods. In the special case of a quartic 
potential with the "wrong" sign, we use functional integration techniques to prove that 
the PT-symmetric model is equivalent to a hermitian matrix model with an anomaly 
for all values of N, as in the one-component case [9l [10]. Interestingly, the anomaly 
vanishes to leading order in the large- limit. 

Although the construction of PT-symmetric matrix models has proved to be 
relatively straightforward, the construction of PT-symmetric models with fields 
transforming as vectors under 0{N) is more difficult technically. Nevertheless, the 
development of scalar field theory models with vector symmetry is crucial to the possible 
relevance of PT-symmetry in particle physics. Only models with quartic interactions 
have so far proved tractable. This progress on quartic models with components 
is built upon recent work on the relation of the one-component — Ax^ model to its 
equivalent Hermitian form [9l [10], as well as recent work on the relation of 0{N)- 
symmetric Hermitian models to one-component PT-symmetric models [H]. As in the 
single- component and matrix cases, the PT-symmetric model with 0{N) symmetry 
and quartic interaction also proves to have a Hermitian form for all values of A^. The 
Hermitian form of the PT-symmetric 0{N) model allows a technically straightforward 
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construction of the large- iV limit, which can in turn be compared with simpler methods 
that lead to essentially the same result at leading order. The constraint-field method 
is particularly notable, apart from its simplicity and familiarity, because it generalizes 
to provide the form for the effective potential of PT-symmetric scalar field theories in 
the large-A^ limit. This effective potential in turn implies asymptotic freedom in four 
dimensions, a property long suspected to hold in a renormalizable PT-symmetric scalar 
field theory in four dimensions. 

This paper is organized as follows: Section 2 develops the formalism required to 
treat PT-symmetric matrix models, and section 3 analyzes the properties of the ground 
state of such models using WKB methods. Section 4 treats the special case of the TrM^ 
model, which has a simple Hermitian dual representation. In section 5, a large class of 
PT-symmetric models with components are described, including vector models with 
0{N) symmetry. Section 6 shows that models in this class have simple Hermitian dual 
representations. In section 7, the large- limit of the 0{N) model is derived using three 
different methods, while section 8 applies one of these methods to PT-symmetric scalar 
field theories in the large- iV limit. 



2. Formalism for Matrix Models 



The techniques for solving Hermitian matrix models are well-known. The solution for 
all N of the quantum mechanics problem associated with the Euclidean Lagrangian 

2 



L = -Tr 
2 



dM 



+ —TrM^ 
N 



(1) 



where M is an x Hermitian matrix was first given by Brezin et a/ [H]. The ground 
state ip is a, symmetric function of the eigenvalues A-,- of M. The antisymmetric wave 
function defined by 



0(Ai,..,AAr) 



j<k 



Afc) 



V^(Ai,..,AAr) 



satisfies the Schrodinger equation 
1 92 



E 



2dX) 



9 ,4 



(2) 



(3) 



where E^^Hs the ground state energy scaled for the large-A^ limit. This equation 
separates into A^ individual Schrodinger equations, one for each eigenvalue, and the 
antisymmetry of determines N'^E^^^ as the sum of the A^ lowest eigenvalues. 

Here we solve the corresponding problem where the potential term is PT-symmetric 
but not Hermitian. As shown by Bender and Boettcher [1], the one- variable problem 
may be solved by extending the coordinate variable into the complex plane. This implies 
that for PT-symmetric matrix problems, we must analytically continue the eigenvalues 
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of M into the complex plane, and in general M will be normal rather than Hermitian. 
We consider the Euclidean Lagrangian 



L 




J^Trmr (4) 



with g > 0. Making the substitution M UAU^, with U unitary and A diagonal, we 
can write L as 

2 




1 



j,k ^ 



\2 




where 



— = -tU^—. (6) 
dt dt ^ ' 

In the analysis of conventional matrix models by Brezin et ai, a variational argument 

shows that the ground state is a singlet, with no dependence on U. Because the Aj's are 

in general complex for PT-symmetric theories, this argument does not apply. However, 

in two cases we can prove that the ground state is indeed a singlet: for p = 2, which 

is trivial, and for p = 4, where the explicit equivalence with a hermitian matrix model 

proven below can be used. Henceforth, we will assume that the ground state is a singlet, 

but our results will apply in any case to the lowest-energy singlet state. 

We have now reduced the problem of finding the ground state to the problem of 

solving for the first states of the single-variable Hamiltonian 

= (7) 

This Hamiltonian is PT-symmetric but in general not Hermitian. The case p = 2 is the 
simple harmonic oscillator. For p > 2, the Schrodinger equation associated with each 
eigenvalue may be continued into the complex plane as explained in [1] . We exclude the 
case p < 2, where PT symmetry is spontaneously broken and the eigenvalues of H are 
no longer real. 



3. Ground State Properties of Matrix Models 

As with Hermitian matrix models, the ground state energy is the sum of the first A^ 
eigenenergies of the Hamiltonian H. In the large A^ limit, this sum may be calculated 
using WKB. A novelty of WKB for PT-symmetric models is the extension of classical 
paths into the complex plane. This topic has been treated extensively for single- 
component models f^ \n\. 

We define the Fermi energy Ep as the energy of the A^'th state 

N=^ [ dpdX 9 [Ep - H{p, A)] (8) 

271 J 
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where the path of integration must be a closed, classical path in the complex j9 — A 
plane. In order to construct the large- limit, we perform the rescaling p \fNp and 
A \fN\ yielding 

1 



Hsc{p, A) =^,/-g {^\r 



(9) 



where the scaled Hamiltonian Hsc is related to H by H = NHsc- We introduce a rescaled 
Fermi energy ep given by Ep = Nep, which is implicitly defined by 

1 = ^/ dpdXe [ep-H,,{p,X)]. 

After carrying out the integration over p, we have 

1 = ^ y d\.j2ep + 2g {iXfO [ep + g [iXf] 

where the contour of integration is taken along a path between the turning points which 
are the analytic continuation of the turning points at p = 2. This equation determines 
ej? as a function of g. 

We define a scaled ground state energy E^'^'> by 

1 N-l 



(10) 



111 



E 



(0) 
N 



HE,. 



k=0 



The WKB result for the sum of the energies less than Ep can be written as 
J2Ek = — I dpdX H,,{p, X)e [ep - H^P, A)] 



k=0 



27r 



(12) 



(13) 



dpdX Hsc{p, X)9[ep- Hsdp, A)] 



so that in the large- limit E^^ is given by 

= — 

°" 27r 

The integration over p is facihtated by using equation (fTOj) to insert a factor of ep, giving 

=^F-^ I dpdX [ep - Hsc{p, A)] 9 [ep - H,,{p, A)] . 
The integral over p then yields 

E^^^=ep-^ I dX [2ep + 2g {iXff^ d[ep + g [iXf] . 
on J 

The turning points in the complex A plane are 



(14) 
ving 
(15) 

(16) 



££ I gm(3/2-l/p) 

9 



i/p 



(17) 



A_ 



. 9 . 



i/p 



-in(l/2-l/p) 



We integrate A along a two-segment, straight-line path connecting the two turning points 
via the origin llj. Solving equation (ITO!) for ep, we find 



ep 



IT 



' r(3/2 + i/p) ^ 

sin(7r/p)r(l + 1/p) 



2p 



9 



p+2 



(19) 
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Table 1. The scaled ground state energy E)^ at g = 1 for p = 3 and p = A. The 
finite N resuhs obtained numericaUy rapidly approach the N —t oo limit obtained from 
WKB. 



N 


p=3 


p=4 


1 


0.762852 


0.930546 


2 


0.756058 


0.935067 


3 


0.75486 


0.935846 


4 


0.754443 


0.936115 


5 


0.754251 


0.936239 


6 


0.754147 


0.936306 


7 


0.754084 


0.936347 


8 


0.754043 


0.936372 


oo 


0.753991 


0.936458 



and solving ( 1161) for the scaled ground state energy we have 



p + 2 
3p + 2 



p + 2 
3p + 2 




r(3/2 + i/p) ^ 

sin (vr/p) r(l + 1/p) ^ 



2p 



9 



p+2 



(20) 



For p = 2, this evaluates to -E''^"^ = y (7/2 , in agreement with the explicit result for the 
harmonic oscillator. 

It is very interesting to compare the large- result with results for finite N. The 
low-lying eigenvalues for the Hamiltonian p'^ — {ixY have been calculated by Bender and 
Boettcher in [1] for the cases p = 3 and p = 4; the case p = 2 is trivial. We can use their 
results by noting that the eigenvalues of our Hamiltonian H are related to theirs by 

,2/(p+2) 



9 



-E. 



BB 



2p/(p+2) jY{p-2)/(p+2) • ^"^^^ 

Results for p = 3 and 4 and small values of are compared with the large- limit in 
Table 1. The energies for finite values of rapidly approach the N ^ 00 limit. The 
approach to the limit appears monotonic in both cases, but with opposite sign. 



The expected value of (TrM) for large A^ is given by 



N-l 



{TrM) = ^ (A,) = — / dpdX Xe [Ef - H{p, A)] . 
j=o 27r J 



(22) 



Calculations of higher moments {TrM'"') are carried out in the same manner. Upon 
rescaling, we find that {TrM) grows as N^/'^^ and the scaled expectation value is given 
by 



lim 



1 



AT^oo A^3/2 



(TrM) 



1 

2^ 



dpdwe [tF-Hsc{pA)] 



(23) 
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which reduces to 



^ = }- j dX \^2eF + 2g {iXyO [2eF + 2g {iXf] . (24) 
Using the same two-segment straight hne path as before, we find that 



/ vr \ cos (A 



(^) 



2 

p+2 



r(3/2 + i/p) 
r(i + i/p) 



p+2 



r (1 + 2/p) 

r(3/2 + 2/p)- 



(25) 



For p = 2, /X = 0, as expected for a harmonic oscillator. For p > 2, the expectation value 
fi is imaginary because (A-,) for each eigenstate of the reduced problem is imaginary jl]. 
For p = 3, fi = — 0.52006i. For p = 4, fi = — 0.772539z. In the limit p ^ oo, fi goes to -i. 
This behavior is easy to understand, because in this limit, the turning points become 
degenerate at —i. 

4. Special case of TrM^ 

For the case of a TrM^ interaction, we can explicitly exhibit the equivalence of the PT- 
symmetric matrix model with a conventional Hermitian quantum mechanical system. 
As in the single- variable case, there is a parity-violating anomaly, in the form of an extra 
term in the Hermitian form of the Hamiltonian, proportional to h. We show below that 
the anomaly term does not contribute at leading order in the large- iV limit. 

The derivation of the equivalence closely follows the path integral derivation for the 
single- variable case [9|, [10] . The Euclidean Lagrangian is 

L = -Tr( ^] ' + -m^Tr - ^TrM^ (26) 
2 \ dt J 2 N ^ ^ 

and the path integral expression for the partition function is 

Z = J [dM] expi^- J dt l| . (27) 
Motivated by the case of a single variable, we make the substitution 



M = -2iVl + iH (28) 

where H is an Hermitian matrix. Because M and H are simultaneously diagonalizable, 
this transformation is tantamount to the relation 



A^- = -2i^l + ihj (29) 

between the eigenvalues of M and the eigenvalues hj of H. The change of variables 
induces a measure factor 

[dM] = (30) 

where the functional determinant depends only on the eigenvalues of H . The Lagrangian 
becomes 

2 1 + iH ^ ' N ^ ' ^ ' 
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at the classical level. However, following [TU], we note that in the matrix case the change 
of variables introduces an extra term in the potential of the form 

2 



' d 


(dhA 


dhj 


[dXj_ 



which can be written as 

32 ^ 1 



AV 



Tr 

32 



1 + iH 



The partition function is now 

[dH] 



det 



VTTTh 



exp 



dtL 



(32) 



(33) 



(34) 



where 



1 {dH/dty , 

-I r 2m I r{l + iH) 

2 1 + iH V ^ ; 



16g 



Tr{l + iHf 



32 VI 



N ^ '32 W + iH^ 

We introduce a hermitian matrix-valued field 11 using the identity 



det 



J [dU] expl^- J dtTr 



- (1 + zH) (n - Ho) 



(35) 
(36) 



where Ho = (iH + 1/4^ / (1 + iH). Dropping and adding appropriate total derivatives 
and integrating by parts yields 



[dH] [dU] exp 



- / dtL' 



(37) 



where 



L' = Tr 

+ Tr 



2m2(l + iH) - 16-^(1 + iH)^ + ^ + 



U{l + zH)--U 



The integration over H is Gaussian, and gives 

N 



[dH] exp < — / dtTr 



Q4g 



4 



After the rescaling 11 J32g/NIl we have finally 



[dU] exp < — / dtTr 



1 



-n 

4 



(3^ 



(39) 



This proves the equivalence of the PT-symmetric matrix model defined by 

2 



L = iTr 

2 \ dt 



+ -m^TrM^ - —TrM^ 
2 N 



to the conventional quantum mechanics matrix model given by 

2 



L' = -Tr 
2 



diV 
dt 



TrU - m'^TrU^ + — TrH^ 
A^ 



(40) 



(41) 



(42) 
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This equivalence implies that the energy eigenvalues of the corresponding Hamiltonians 
are the same. This could also be proven using the single-variable equivalence for the 
special case of singlet states, but the functional integral proof encompasses both singlet 
and non-singlet states at once. The equivalence of these two models also allows for an 
easy proof of the singlet nature of the ground state. Standard variational arguments 
show that the ground state of the Hermitian form is a singlet. The direct quantum 
mechanical equivalence of the single-variable case is then sufficient to prove that the 
ground state of the PT-symmetric form is also a singlet. 

As in the single-variable case, there is a linear term of order h appearing in the 
Lagrangian and Hamiltonian of the Hermitian form of the model. This term represents 
a quantum mechanical anomaly special to the TrM^ model. To determine the fate of 
the anomaly in the large- limit, we construct the scaled Hamiltonian of the Hermitian 
form in exactly the same way as for the PT-symmetric form. It is given by 




indicating that the effect of the anomaly is absent in leading order of the large-iV 
expansion. One easily checks for the m = case that the Hermitian form without the 
linear term reproduces the PT-symmetric prediction for E'^ at p = 4. 



5. 0{N) Vector Models 



The analysis of the 0(A^)-invariant PT-symmetric model with a quartic interaction is 
similar to that of the TrM^ matrix model. Consider a model with Euclidean Lagrangian 
given by 

2 



N r 



E 



1 



H — m X,: — \x, 
2 




(44) 



where g and A are non-negative. When (7 = 0, we have decoupled one-dimensional 
systems; for A = 0, we have a model with 0{N) symmetry. When both g and A are 
non-zero, the model has only an Sn permutation symmetry. From the standpoint of 
PT symmetry, the interaction terms can be considered as members of a family of PT- 
invariant interactions 



N 



\2p 




(45) 



which are invariant under PT symmetry. This class of models is well-defined for 
p = q = 1, and must be defined for p, g > 1 by an appropriate analytic continuation of 
the Xj as necessary [1]. 

It is convenient to consider this model as a subset of a larger class of models, with 
a Lagrangian of the form 



E 



[dtx 



"2 , 2 2 

3J + 



j,k=l 



(46) 
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The classical stability of the potential for large Xj is governed by the eigenvalues of A. 
For the model of particular interest to us, 

A = XI + gP (47) 

where P is the one-dimensional projector 



1 

N 



( 



V 



1 \ 



(48) 



satisfying P^ = P. The decomposition A = X {I — P) + {g + X) P shows that A has one 
eigenvalue g + X and — 1 eigenvalues with value A. The eigenvalue (7 + A is associated 
with variations in x^, i.e., variations in the radial direction. 



6. Equivalence of PT-symmetric vector models to Hermitian models 



We will analyze the case where all eigenvalues of A are positive using functional 
integration. With the substitution 

—2i^Cj + iipj 



(49) 



familiar from the one-component case, Le becomes 

"1 (dt^P.f 



E 



2m {cj + ii/jj) 



y [2 (cj+iijj) 
- 16^ Ajfc(cj + ii'j){ck + iipk)- 

jk 

The generating function for the model is given by 



n 



det(co + iipj 



exp 



dt 



t32 



1 



(50) 



(51) 



where the change of variables has generated both a functional determinant and 
additional term, formally of order in the action. As pointed out in [10], both terms 
are required to obtain correct results in the functional integral formalism. 
The functional determinant may be written as 
1 



n 



det [-y/cj + itpj 



J Y[ [dhj] exp \~J'^^ 



-{cj + iijj) [hj - 



iijj + 1/4' 



(52) 



which introduces a new set of fields hj. The derivation proceeds as in the single- 
variable case. After integration by parts on the hjtpj terms, and adding and subtracting 
total derivatives, the functional integral over the ipj fields can be carried out exactly. 
The integral is both local and quadratic, and requires that the matrix A have positive 
eigenvalues for convergence. The result of this integration is 



exp 



dtL 



H 



(53) 
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where Lh is given by 

1 

64 



-/i? + L - 2m^ 
2 J -I 



-hl + hk- 2m^ 



(54) 



After discarding total derivatives, we obtain 



h^h 



(55) 



which gives the Hermitian form for our general PT-symmetric model with N fields. 
In the particular case we are interested in, we have 



A 



-1 



' g+X 



(56) 



The Lagrangian may be written as 



1 



64A 



E 



9 



Am' 



64A^A (g + X) 



(57) 



It is helpful to immediately rescale all the fields as hj 



■IE 



9 



VMXhj-. 
-V2Xj2hj 



N(g + X) 




(58) 



At this point, the Sn permutation symmetry is still manifest, and it clear that the field 
J2j hj plays a special role. 

In order to understand the strategy for rewriting the model in a form in which the 
hmit A — > can easily be taken, it is useful to work out exphcitly the case oi N — 2 
first. It is apparent that a rotation of the fields will be desirable. We define suggestively 
new fields a and tt given by 

/l2 = ^ ((7 - tt) ■ 

After some algebra and the rescaling 



(59) 



19 + ^ 



we arrive at 



^E 



X 



2 2 



(60) 



Am^ 
^A' 



-TT 



2(^ + A)a^ + ^7r^ 



^ + A 



+ (8^ + 12A) crV^ - 2^g + Xa. 



(61) 
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Notice the natiiral hierarchy between the masses for \ <^ g. The 0(2) symmetric hmit 
of the original PT-symmetric model is obtained in the limit A ^ 0, where we have 

Le = + ^TT^ - mV^ + 2^(T^ + Sga^n'^ - 2^a. (62) 

The field tt has no mass term, indicating its relation to the angular degrees of freedom in 
the original Lagrangian. However, radiative corrections generate a mass for the tt field 
via the the (T^tt^ interaction. As in the one-component case, there is a linear anomaly 
term, but only for a. 

We now turn to the more difficult case of the A — > hmit for arbitrary N. As 
before, we introduce a field a defined by 

cr = -]= hj (63) 



N 

as well as a set of iV — 1 fields tt^ with k — 1,..,N — 1 related to the hj fields by a 
rotation so that a"^ + n'^ — . Each field hj can be written as 

hi = ^L(7 + hi (64) 



where J2j hj = 0. This property is crucial in eliminating a term in which diverges 
as A~^/^ as A ^ 0. The Lagrangian now can be written as 




{a' + 7f2) + 4A ^ /ij 



ox 2 



-AJ [a^ + Tf^-^l -V2XNa. (65) 



The rescaling a '^{g + X) / Act plus some careful algebra yields the A — > limit as 

^E^-o- +-7T -ma + —a + tt - ^2gNa (66) 

which agrees with our previous result for = 2 , and agrees with the known result 
for a single degree of freedom if we take A^ = 1 and drop the vf field altogether. This 
is a Hermitian form of the PT-symmetric anharmonic oscillator with 0{N) symmetry, 
derived as the limit of a PT-symmetric model with symmetry. The Hermitian form 
has several novel features. Note that both the Sn and 0{N) symmetries are no longer 
manifest, but there is an explicit 0(A^ — 1) symmetry associated with rotations of the tt 
field. As in the N — 2 case, there is no mass term for the tt field. Furthermore, there is 
no (tt^)^ term, although there is a tt^o"^ interaction. The anomaly term again involves 
only (7, and breaks the symmetry a —a possessed by the rest of the Lagrangian. 
Analyzing the Lagrangian at the classical level, we see that if > , the a field is 
moving in a double-well potential, perturbed by the anomaly so that {a) > 0. On the 
other hand, if < 0, a moves in a single- well anharmonic oscillator, again with the 
linear anomaly term making (a) > 0. In either case, the vr^o"^ interaction will generate 
a mass for the tt field. All of this is consistent with the association of a and tt with 
the radial and angular degrees of freedom, respectively, in the original PT-symmetric 
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model. This equivalence between PT-symmetric and Hermitian forms may be compared 
to the results of [11], where a somewhat different equivalence is derived. In that work, 
the generating function for a Hermitian theory is shown to be equivalent to a sum 
over generating functions for a class of single- component PT-symmetric models, with 
each element of the class representing a different angular momentum. As we discuss 
below in the context of the large- limit, both approaches lead to an anomaly term 
with a linear dependence on the angular momentum quantum number /. 



7. Large- Limit of Vector Models 

We will defer a more detailed discussion of this model for finite A^, and turn to its 
large- A^ limit. One more rescaling a \fNa gives the Lagrangian 

Le = y cr^ + - Nm^a^ + AgNa^ + Wga'^n^ - N^ga. (67) 

We see that the anomaly term survives in the large- A^ limit, unlike the matrix model 
case |1] . After integrating over the A^ — 1 components of the tt field, we have the large- A^ 
effective potential K// for a: 



V,ff/N = -m'a' + Aga'' + -^?>2ga^ - ^2ga. (68) 

It is striking that the anomaly term has virtually the same form as the zero-point energy 
of the vf field. The anomaly term breaks the discrete a —a symmetry of the other 
terms of the Lagrangian, and always favors cr > 0. The effective potential has a global 
minimum with o positive for w? > 3 2^^^g'^^^. For < 3 2^^^g'^^^, there does not appear 
to be a stable solution with a > 0, and a = is the stable solution to leading order 
in the 1/A^ expansion. This change in the behavior of the effective potential as is 
varied is not seen in the corresponding Hermitian model [H], and indicates a need for 
care in analyzing the model. Based on our preliminary analysis of the Hermitian form 
for finite A^, we believe that this behavior is associated with the large- A^ limit, and does 
not indicate a fundamental restriction on m^. 

The large-A^ effective potential was derived from a Lagrangian with unusual 
properties, associated with the Hermitian form of the original model. It is therefore 
surprising that, once the form of the large- A^ effective potential is known, it can be 
derived heuristically in a more conventional way. We start from the 0(A^)-symmetric 
Lagrangian 



^ rl 1 
and add a quadratic term in a constraint field p 




(69) 
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1 _ .2 , , 2I , 4A^P^ Nrn^p Nm 



9 9 1% 



yielding 

N 

i=i 

If we integrate over Xj in a completely conventional way, we obtain the large- iV effective 
potential 

K„/iv = l^-!!!V + y^+!^l. (72) 

9 9^ 16c/ 
This is essentially identical to our previous expression after identifying p = ga^. 
However, we lack a fundmental justification for this approach. We know that great 
care must be taken in specifying the contour of integration in typical PT-symmetric 
models, yet the xj fields were integrated over quite conventionally. If this approach has 
validity, it seems likely that the choice of integration contours for p and x is crucial. 
However, only the saddle point matters to leading order in so it is possible for this 
heuristic derivation to be correct even though we lack a direct, complete treatment of 
the original PT-symmetric model. 

There is another approach to the effective potential that sheds some light on the 
role of angular momentum in PT-symmetric vector models. Let us take as our starting 
point the Hamiltonian for the PT-symmetric vector model after the introduction of the 
constraint field: 



N 



+ (73) 

9 9 



where for simplicity we have dropped the constant term. The reduced Hamiltonian 
for the radial degree of freedom can be written as 

(N + 2l-l)(N + 21-3) , 2 4Arp2 Nm^p 
2 or'^ Sr"' g g 

where / is a non- negative integer [13]. Rescaling the radial coordinate r — N^/'^r leads 
to a potential proportional to N which is a function of both r and p and a kinetic term 
which is of order It is easy to minimize the potential as a function of r; the final 

result, after setting Z = 0, is identical to the expression for Vejf/N as a function of p in 
the large- iV limit. Thus we see that this radial formalism yields results for the ground 
state energy equivalent to other approaches at leading order in the large- iV expansion. 
Alternatively, one can take the angular momentum quantum number / to be of order 
A^. Elimination of the r variable then leads to an effective potential of the form 

ANp^ Nrri^p 



K// = (A^ + 2/) J2p + ^ (75) 

^ 9 9 



which displays the /-dependent anomaly term first observed in [TT]. It is also possible to 
show the equivalence of the radial formalism directly, without introducing the composite 
field p. Note that the radial approach demonstrates that the angular momentum term 
makes a positive contribution to the ground state energy in the PT-symmetric case, 
exactly as it does in the Hermitian case. 
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8. PT-symmetric field theory 

If we boldly apply the constraint field approach to a PT-symmetric field theory with a 



Models of this type were rejected decades ago [14j because of stability concerns at both 
the classical and quantum levels, although there were early indications that such theories 
were in fact sensible [I5]. Within the framework of PT-symmetric models, such stability 
issues cannot be addressed without a detailed understanding of the contours used in 
functional integration. However, it is straightforward to check that renormalization of 
g m. d = A gives an asymptotically free theory, with beta function (3 = —g^j^'n'^ in the 
large- limit. If PT-symmetric scalar field theories exist in four dimensions and are 
indeed asymptotically free, the possible implications for particle physics are large, and 
provide ample justification for further work. 

The authors would hke to thank Carl M. Bender for many useful discussions; we also thank Stefan 
Boettcher for providing the numerical data used in constructing Table 1. MCO gratefully acknowledges 
the support of the U.S. Department of Energy. 
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